function model = kdaclassify(stego, cover, ktype)
    if nargin < 3
        ktype = 'Gaussian';
    end
    M1 = size(stego, 1);
    M2 = size(cover, 1);
    fea = [stego;cover];
    gnd = [ones(M1, 1);zeros(M2, 1)];
    options.KernelType = ktype;
    options.t = 1;
    [eigvector, eigvalue, K] = KDA(options, gnd, fea);
    model = struct('eigvector', eigvector,'eigvalue',eigvalue,'options',options,'kernel',K);
end
